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Abstract 



?7c(2980) production in 77 interactions has been detected via its decays into 
K^K^n^, K+K-K+K- and K+K~7T+7r- in the data taken with the DELPHI 
detector at LEPl and LEP2 energies. The two-photon radiative width averaged 
over all observed decay channels is P^^ = 13.9 ± 2.0 (stat.) =b 1.4(syst.) ± 2.7 
(BR) keV. No direct decay channel rjc — > 7r"'"7r~7r"'"7r~ has been observed. An 
upper limit P^^ < 5.5 keV at 95% confidence level has been evaluated for this 
decay mode. 
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1 Introduction 

Among 77 induced final states, those with exclusive meson resonance production play 
an important role, since the measurement of the production cross-section and the corre- 
sponding radiative width provide information on the quark-gluon structure of the investi- 
gated particle. Among these final states, those with mesons built up of heavy quarks are 
particularly interesting since such mesons can be described with nonrelativistic models. 
In particular, a precise measurement of the two-photon partial width for charmonium 
states would provide valuable information on QCD corrections to cc quarkonium. 

The very first estimations of the rjc partial width, r^^(?7c), were obtained from its 
ratios to the known widths for ^ f^^fJ^~ and rjc —>■ gg giving values of 8 keV and 4 
keV respectively [1]. Different models and corrections were applied to them later, giving 
values from 3 to 14 keV, see [2] and references therein. An even bigger discrepancy is 
observed between values obtained by numerous experimental groups. Among them there 
are many experiments where two interacting photons radiated by electron and positron 
beams couple to this resonant state, [3]- [12]. The results for r^^{r]c) range from 4 keV 
to 27 keV. 

In this paper we report on the production and decays of the rjc resonance using data 
collected by the DELPHI detector during the period 1994-1999 corresponding to a range 
of centre-of-mass energies from 90 GeV up to 202 GeV and an integrated luminosity of 
jC — 531 pb~^. The aim of this analysis was to determine the radiative width of the rjc 
resonance separately for each decay channel, using the production process: 

e+e" ^ e+e"?7c(2980) (1) 

on four-body final states where a distinct signal of the rjc resonance has been observed. 
To increase the sensitivity for rjc production, we do not require information on the polar 
angle of the scattered electrons (no tag mode). The superiority of LEP with respect to 
previous experiments is the higher energy and resulting higher production cross-section 
for this reaction. 

We have analysed the following exclusive final states: 



rjc ^ X°X^7r^ (2) 

77c ^ K+K-K+K- (3) 

rjc K+K-n+n- (4) 

77c — > 7r''"7r~7r''"7r~ (5) 



2 Detector 

A general description of the DELPHI detector can be found elsewhere [13]. The 
main features relevant to this analysis are particle tracking and identification. Due to 
the low momenta of the decay products, their identification is based on measurement of 
ionization losses (dE/dx) in the Time Projection Chamber (TPC). The particle momenta 
are determined from track reconstruction and make use of the Vertex Detector, the Inner 
and Outer Detectors and the TPC. The tracks with lower polar angles are reconstructed 
in Forward Chambers A/B. 

The single track trigger efficiency, expressed in terms of transverse track momentum, 
has an influence on the overall efficiency of flnal states produced in 77 collisions where 
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the hadrons have rather low momenta. Having four particles in the final state, originating 
from the decay of a relatively heavy r^c (2980) resonance, results in a large trigger efficiency 
for an event according to the formula: 

£ev = l-{l- ei) X (1 - £2) X (1 - es) X (1 - e4) (6) 

where E^v is the total trigger efficiency for an event and the is the single track efficiency 
which depends on the transverse momentum. A brief description of the trigger system is 
presented in [14,15]. 



3 General Data Selection 

Data were taken only from running periods when the TPC was fully operational thus 
ensuring good particle identification. There was no requirement on detecting either scat- 
tered electron. Candidates for the ?7c(2980) decay channels (2)-(5) were selected by re- 
quiring: 

• exactly four charged particle tracks with zero total charge, coming from the primary 
interaction region or two tracks originating from the primary vertex and two tracks 

originating from a secondary vertex, 

• the track impact parameters measured with respect to the z-axis (beam axis) to be 
smaller than 10 cm and those measured in the plane perpendicular to the z-axis 
smaller than 4 cm, 

• the momentum of each particle to be larger than 0.1 GeV/c, 

• the square of the total transverse momentum, (EpJ)^, of charged particles to be less 
than 1.0 (GeV/c)2, 

• each track to pass through the TPC, 

• the total detected energy of charged particles to be less than 10 GeV, 

• no particles identified as electrons or muons by the standard lepton identification 

algorithms, 

• the track lengths to be longer than 30 cm, 

• the total energy deposit in the electromagnetic calorimeter from neutral particles to 
be less than 3 GeV, 

• the charged particles to have polar angles between 20° and 160°. 

Additional criteria which are specific to particular channels are discussed in the next 
section. 

All experimental requirements used in the analyses presented below were chosen to be 
the same for all data sets corresponding to various beam energies. 

4 Analysis 

In 77 events almost all the available energy and momentum is carried away by the 
electron and positron which are scattered at very small angles. Therefore the (Sj9^)^ 
distribution of the hadronic system is peaked at low values, as shown in Fig.l. To 
suppress background events which do not originate from 77 collisions, the total transverse 
momentum squared of hadrons in the exclusive process (1) should be smaller than 0.04 
(GeV/c)2. 
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In order to calculate the acceptance and detection efficiency, a Monte Carlo generation 
program has been used, with the full kinematics of a system produced in 77 interactions. 
All kinematical variables necessary for the description of the two-photon processes were 
generated using algorithms taken from the package described in [16]. The matrix ele- 
ment, factorized into the flux of quasi-real transverse photons and a covariant amplitude 
describing both the two-photon rjc production and its decay, has been implemented [17]. 
For a better understanding of the rjc four pion decay mode we have also determined the 
efficiency for rjc p^p^ — > 7r"'"7r~7r"'"7r~ with a specific symmetrized matrix clement [17]. 
The Monte Carlo generated events were passed through the standard DELPHI detector 
simulation procedure [13]. 

An additional factor contributing to the overall efficiency comes from the trigger ac- 
ceptance. The trigger simulation following the cuts used for r]c selection in the real data 
has been applied to events after detector simulation. An event was accepted according to 
a weight calculated on the basis of the single track efficiency, parameterized as a function 
of the transverse momentum, pti and ranges from 20% for pt= 0.5 GeV/c to about 95% 
at pt— 2 GeV/c [14,15]. Owing to the relatively large mass of the ric resonant state, the 
overall trigger efficiency per event was about 90% for channels with pions and about 85% 
for the K+R-K+K- final state. 

The total efficiency was calculated bin-by-bin in invariant mass by comparing the 
generated invariant mass distribution with that obtained from the detector simulation 
after the selection cuts and trigger acceptance. The efficiency for each decay mode as 
a function of the invariant mass is shown in Fig.2. It should be noted that particle 
identification was essential for all the channels analysed and was based on dE/dx energy 
loss measurements [13]. 

4.1 r/e ^ K^K^Tr^ 

For the decay chain r]c — >■ K^K^n^ ti^ti~ K^n^ the K'^ tt+tt" decay is identified 
by taking advantage of the relatively large decay length (cr = 2.68 cm). Therefore, 
candidates for this decay mode had to have one secondary vertex reconstructed using 
an algorithm which takes pairs of oppositely charged particle tracks, intersecting them 
and determining a secondary vertex. Both momenta arc recalculated with respect to the 
new decay vertex and an invariant mass is computed. The resulting candidate mass 
distribution is shown in Fig. 3, where clear evidence of a signal is seen. Only events 
with an invariant mass of the two pion candidates, originating from the secondary vertex, 
in the range from 0.45 GcV/c^ to 0.55 GeV/c^ have been taken for further analysis. Of 
the other two particles which originate from primary interaction region, one is identified 
as a kaon in 80% of the events selected with one secondary vertex. Hence the crucial 
criterion for this decay final state selection is the reconstruction of the decay vertex. 

4.2 r]c K+K-K+K- 

Additional requirements for this decay channel are that at least three particles must 
be identified as charged kaons and there arc no secondary vertices. Only kaons with the 
probability of identification greater then 0.5 were considered. The dE/dx distribution for 
all identified particles after the general data selection is plotted in Fig. 4 with an insert 
for the distribution of those originating from r/c (2850 MeV/c^ < M{K+K-K+K~) < 
3150 MeV/c^). For events from this rjc mass region, points originating from the rising 
part of this distribution unquestionably correspond to kaons whereas the horizontal part 
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may also contain pions from background events and kaons from signal. 
A scatter plot (not shown) of the invariant mass of K'^K~ combinations does not indicate 
any intermediate (fxf) state. Prom a fit to the invariant mass distribution, the number of 
signal events is estimated to be about 46. 

Since the average particle momentum is particularly low in this channel, a strong ef- 
fect could be expected in the invariant mass spectrum resulting from the single track 
efficiency of the trigger that might produce a fake signal due to the small efficiency at 
threshold. This has been checked on rjc — > K~^K~K~^K~ events which were generated 
according to the 77 flux (no resonance shape has been assumed) and then decayed ac- 
cording to phase-space. These events were then passed through the trigger and detector 
simulations. No signal resulting from the trigger activity on the low mass side nor from 
the experimental cuts on the other was observed in a region of the invariant masses 
around 3 GeV/c^, corresponding to the rjc signal. The relatively low background at the 
K^K~K^K~ invariant mass threshold is explained mostly by the low acceptance and 
less by the decreased trigger efficiency. The trigger efficiency, as described in previous 
section turned out to be around 85% at a mass of 3 GeV/c^. 



4.3 r]c K+K-ir+ir- 

Given the branching ratio for rjc decay into K'^K~7r'^7r~ , BR=2.0±0.7%, and the 

detector efficiency determined using criteria presented below (Fig. 2), a significant signal 
(of about 4 events per 1 keV of rjc radiative width) would be expected in this channel. In 
order to select these events it was required that one of the particles was identified as a kaon 
with probability >0.5 and two of the three remaining particles should satisfy selection 
criteria for pion identification with probability >0.5. The identification was based on 
dE/dx energy losses. All events corresponding to the K^K^n^ signal, described in section 
4.1 have been subtracted from the selected sample. Since the data sample obtained 
may still contain K'^K^t^^ events with no reconstucted secondary vertex, the invariant 
mass Mjfe of the two opposite sign particle combinations (excluding the identified kaon) 
was calculated and events removed if one of the two Mj^ masses satisfied the condition 
\Mkq — Mik\ < 50 MeV/c^. From the data collected by the DELPHI detector during the 
period mentioned in the first section a signal of about 42 events is obtained. 

The intermediate states of rjc — > K'^K~7r'^7r~ decay via one or two K*^{S92) have not 
been observed. 

4.4 ijc — ^ 7r'^7r~7r''"7r~ 

The observation of an rjc — > 7r''"7r~7r''"7r~ decay mode reported by numerous exper- 
imental groups remains controversial. This decay has been found by MARK III [4], 
DM2 [9], TASSO [7] (where the last one did not distinguish between the An and the p^p^ 
decay channels). Among more recent experiments this final state has been observed in 
EES [18]. None of the LEP experiments confirm this decay mode providing only an upper 
limit [11]. Good particle identification is very important since the 7r"^7r'"7r"^7r~ final state 
can be confused with the K'^K~7r^7i^ decay. 

In addition to the general selection and the stringent cut on the total transverse 
momentum, it was also required that all particles were pions with the single track iden- 
tification probability > 0.5, that only one well reconstructed vertex was found and that 
each track had to have at least one hit in the Vertex Detector. The final selected sam- 
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pie consists of ~3600 events and shows no enhancement around the nominal mass of 
the r)c resonance, see Fig. 5a. Using the PDG values [19] for the rjc parameters, more 
than 60 events would be expected in this channel. An upper limit of 26 events at 95% 
confidence limit has been calculated. The above, standard selection criteria lead to an 
invariant 7r+7r^7r+7r~ mass distribution with a large background that may shadow the 
signal. Further tightening of the total transverse momentum squared cut from 0.04 to 
0.004 (GeV/c)^ and the identification probabihty from 0.5 to 0.8 reduces the number of 
observed events to about 600 but still the invariant mass distribution shows no evidence 
for rjc, see Fig. 5b. To avoid a selection bias resulting from the low efficiency for the iden- 
tification of four pions another selection was performed in which only three particles were 
identified as pions with probability > 0.5, leaving remaining cuts like in the standard 
selection, again resulting in no enhancement at r)c invariant mass region, see Fig. 5c. A 
search for the intermediate decay mode, through an analysis of the two-dimensional 
plot of the invariant mass of one tt+tt" system versus that of the remaining vr+vr" pair 
has been also performed. Events from the mass window were selected and used for 
the calculation of the 7r'*"7r~7r''"7r~ invariant mass spectrum. Since the rjc signal was not 
seen, these events were attributed to non-resonant p^p*^ vector mesons production. 



5 Results 

Experimentally one measures directly the invariant mass {W^^) distribution of the 77 
system, 

where f denotes one of the investigated decay modes. Given the detector efficiency Sf, 
the integrated luminosity C, and flux L^^ of the two interacting photons parametrized 
by well known equivalent photon approximation formula, the invariant mass distribution 
can be converted into two-photon cross section multiplied by corresponding branching 
ratio: 

The efficiency was calculated dividing bin-by-bin the simulated invariant mass distri- 
bution for events that passed all the cuts in the mass interval 2.5-4.0 GeV/c^ by the 
invariant mass distribution for the generated events. It should be noticed that both the 
77 flux and the invariant mass efficiency distribution modify the background-to-signal 
ratio measured in the side-bands of the AN/ AW^^^ distribution. 

In order to determine the value of the rjc radiative width, the 77 invariant mass cross- 
section has been fitted to the Breit Wigner distribution of the form 

BW{r,„ r,„„ W,,) = 87r(2J + 1) (^^72^3^yTMfTi; 

describing the rjc production cross-section convoluted with a Gaussian mass resolution 
G{Wy^, cr) together with a background parametrization expressed in terms of polynomial 
function of the third order Ps{Wjj): 

a^-y^r,Mii) = [BW{T^^, M^,, Ttot, W-y^) + PsiW^^)] ® G{W^^, a) (10) 

According to eq.(8) and eq.(9) the fit determines the product of the radiative width 
and the branching ratio, the mass of the resonance and the experimental mass resolution 
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a. All these fitted parameters have been determined separately for each data sample 
because some of them explicitly depend on the energy (two-photon flux) and others on 
the period of the data collection (efficiency). The total width Vtot of the resonace has 
been fixed to value obtained by other experiments [19]. 

The width of the mass resolution distribution obtained from the above fit coincided 
within ±10% with that obtained from the simulated sample. 

The final a^^ plots are average distributions from different samples. The resulting 
cross-sections multiplied by the corresponding branching ratios for r^^ ~^ K^K^n^, rjc — > 
K^K~K~^K^ and rjc K^-k^tt' are presented in Fig.6. 

A major contribution to the systematic uncertainty originates from the cuts variation 
(about 44% of the total systematic error) and different fit ranges as well as the choice 
of binning (about 28%). Since both the branching ratio and F^^ (see formula above) 
cannot be determined simultaneously, we have used branching ratios and corresponding 
uncertaintes obtained by other experiments, summarized in [19] . This sort of uncertainty 
also contributed to the total systematic error (in the amount of 13%). The remaining 
part of systematic error (15%) was due to the uncertaintes of background shape, trigger 
efficiency and integrated luminosity value. 

In summary, the final values of the r]c radiative width for the three decay channels 
investigated are presented in Table 1: 



final state 


BR (77e final [19]) [%] 




[keV] 




1.5±0.4 


41 


13.3 ± 2.6(stat.) ± 2.0(syst.) ± 3.5(BR) 




2.0±0.7 


42 


14.2 ± 4.9(stat.) ± 2.9(syst.) ± 4.9(BR) 


K+K-R+K- 


2.1±1.2 


46 


16.5 ± 4.3(stat.) ± 2.7(syst.) ± 9.4(BR) 


7r"'"7r~7r"'"7r~ 


1.2±0.4 


< 26 


<5.5 at 95% confidence 



Table 1: The branching ratios taken from PDG, the number of events and radiative 
widths for the particular decay modes. The ^eviVc) is the number of events selected for 
region Af,j^± 150 MeV/c^. In the case of four pions final state only an upper limit has 
been estimated assuming that all events in the rjc mass interval are background. 



For the analysed channels, the results quoted above are the averages of the LEPl and 
LEP2 results. The product F^^- BR in the analysis of the four charged kaon decay channel 
is in agreement, within the large errors, with the result of the ARGUS Collaboration [10], 
which gives 0.231 ± 0.090 (stat.) ± 0.023 (syst.) keV. 

A weighted mean of the radiative width value for the first three channels in Table 1 with 
weights inversely proportional to the total error squared has been determined. 
The result is: 

F.^.^= 13.9 ± 2.0(stat.)±1.4(syst.)±2.7(BR) keV 
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Figure 1: The square of the total transverse momentum of the hadronic system. Points 
represent the real data sample after the general data selection. The histogram shows this 
distribution for dedicated 7]^ production simulation sample with an arbitrary normaliza- 
tion. 
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Figure 2: The average luminosity- weighted efficiences for different rjc decay final states as 
a function of the corresponding invariant mass. In the decay to K^, its branching fraction 
to TT+TT^ has been taken into account. 
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Figure 3: Invariant mass of two particles originating from a secondary vertex (summed 
over all energy samples). 



2 - 



DELPH 


- .'■■■K ' 


3 
2.5 
2 

. ■' P 1-5 
' '. ■ ■ 1 


L ;. T]^ region 


_ • ■'TT ■ ■: , . 


■ ■ ■ • 0.5 


r 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




0.5 1 1.5 ■ 2 


_l 1 1 1 1 1 1 1 1 1 L 


_I 1 1 1 1 1 1 1 1 1 1 L 


_l 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l_ 


0,25 0.5 


0.75 1 1 
P 


.25 1.5 1.75 2 

(GeV/c) 



Figure 4: The dE/dx distribution for particles identified as a pion or a kaon in the data 
events after the general cuts. Most of the remaining tracks consist of protons. The same 
distribution for events from rjc — > K" K~ is shown in the insert. 
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Figure 5: Final invariant mass distributions for the 7r"'"7r t^^ti decay final state. The 
presented distributions are based on: Fig. 5a - the standard selection, Fig.Sb - the stringent 
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